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EFFECTS  OF  PARTIAL  SUPPRESSION  OF  ION  EMISSION 
IN  MODERATE  IMPEDANCE  DIODES 


I.  INTRODUCTION 

The  efficient  generation  of  intense  beams  of  energetic  light  ions  is  a  central  objective  of  the  NRL 
Light  Ion  Fusion  Research  Program.1  Over  the  past  several  years,  focused  ion  current  densities  of  over 
100  kA/cm2  from  terrawatt-level  beams  have  been  achieved  with  magnetically  insulated  radial  diodes  at 
Sandia  National  Laboratories2  as  well  as  with  pinch-reflex  axial  diodes  at  the  Naval  Research  Labora¬ 
tory.3  Using  diodes  of  below  2  ohm  impedance,  ion  beam  efficiencies  of  over  70%  were  achieved  in 
both  of  the  above  configurations4, 5  (i.e.,  over  70%  of  the  power  coupled  to  the  diodes  was  carried  by 
the  light  ions  generated  therein) .  The  question  of  ion  efficiency  is  critical  to  the  goal  of  a  practical  light 
ion  driven  inertial  confinement  fusion  (ICF)  reactor.  As  much  as  possible  of  a  given  pulsed  power 
generator’s  energy  must  be  imparted  to  the  ion  beam  exiting  the  diode  in  order  to  minimize  the 
number  of  beam  sources  necessary  for  successful  pellet  ignition  as  well  as  to  maximize  the  overall  reac¬ 
tor  efficiency.  These  efficiency  considerations  are  well  met  by  the  low  impedance  diodes. 

Unfortunately,  the  high  current  densities  of  the  ion  beams  produced  by  low  impedance  diodes  are 
not  compatible  with  the  focusing  and  transport  systems  presently  under  study  for  bringing  the  beams  to 
bear  on  the  proposed  fusion  targets.6  In  addition,  there  are  strong  arguments  in  favor  of  high  voltage, 
the  use  of  high  impedance  generators  in  present  reactor  scenarios.2  Diodes  matched  to  such  generators 
must  likewise  be  of  high  impedance  to  ensure  efficient  power  transmission  but  such  diodes  are  plagued 
by  relatively  low  ion  production  efficiencies.  These  low  efficiencies  are  a  direct  consequence  of  esta¬ 
blished  diode  theory.  It  has  been  found  semiempirically*  that  the  total  current  flowing  through  a 
pinched-beam  diode  may  be  approximated  by 

/-/,  +  /,  =  9(y2  -  l),/2-| 

where  y  —  1  +  —  q*!^^  .  ^  “  cathode  radius,  D  =  axial  anode-cathode  (A  -K)  gap,  V  —  diode  vol¬ 
tage,  and  m ,  —  ion  (proton)  mass.  Implicit  in  this  formula  is  an  ion-to-electron  current  ratio  given  by 

y-  >  0.5  — y-  (2) 

where  v,  is  the  mean  ion  velocity.9  Thus,  for  a  fixed  voltage,  increasing  the  diode  impedance  translates 
to  decreasing  the  aspect  ratio,  R/D.  That,  in  turn,  results  in  a  decrease  of  the  current  ratio  /,/ /,  and  a 
lowering  of  the  ion  production  efficiency,  IJ(le  +  /*).  This  is  the  crux  of  the  problem  addressed  by  this 
report.  The  essence  of  the  proposed  solution  is  an  attempt  to  decouple  the  impedance  from  its  R/D 
dependence  while  exploiting  the  theorized  efficiency  enhancement  gained  through  a  decreased  gap  size, 
D. 
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For  fixed  values  of  R ,  D ,  and  voltage,  the  net  diode  current  can  be  changed  from  that  predicted  in 
Eq.  (1)  by  restricting  the  emission  of  electrons  from  the  cathode  and  ions  from  the  anode  respectively. 
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As  a  starting  point,  take  the  standard  pinch-reflex  diode  (PRD)  geometry  shown  in  Figure  ).  Typical 
electron  and  ion  trajectories  are  as  sketched.  The  predominant  flow  of  electrons  originates  from  the 
cathode  shank  tip.  Electrons  emitted  there  stream  to  the  thin,  polyethylene  anode  foil  (pinching  radi¬ 
ally  as  they  cross  the  A-K  gap  if  the  critical  diode  current  is  exceeded10).  They  then  pass  through  the 
foil  with  an  appropriate  loss  of  energy.  Finally  they  are  reflected  back  through  the  foil  by  the  large 
azimuthal  magnetic  field,  B9 ,  which  is  present  between  the  foil  and  the  solid  anode  surface.  This 
reflection  (or  "reflexing")  process  continues  as  the  electrons  cascade  toward  the  diode's  central  axis 
where  they  are  funneled  off  through  the  anode  stalk.  The  hydrocarbon  anode  plasma  formed  along  the 
foil  meanwhile  acts  as  a  rich  source  of  protons  which  are  accelerated  in  the  opposite  direction  toward 
the  cathode.  Proton  emission  can  be  expected  everywhere  the  anode  plasma  has  formed  including 
regions  along  the  foil  at  larger  radii  than  the  outer  radius  of  the  cathod  shank.  Ions  emitted  at  radii 
above  the  inner  radius  of  the  cathode  shank  are  useless  for  light  ion  research  for  two  major  reasons. 
First  of  all,  their  trajectories  differ  so  little  from  straight  lines  that,  assuming  a  planar  anode  foil,  they 
will  probably  all  hit  and  be  absorbed  by  the  solid  cathode  shank  and  thus  be  unavailable  for  extraction 
from  the  diode.  Secondly,  their  presence  near  the  cathode  tip  will  further  increase  the  emission  of  elec¬ 
tron  current  with  presumably  detrimental  results  to  the  overall  ion  production  efficiency. 

These  considerations  lead  one  to  conclude  that  significant  benefit  could  be  derived  by  suppressing 
the  emission  of  ions  from  those  outer  regions  of  the  anode  foil  face.  This  is  accomplished  experimen¬ 
tally  by  covering  the  plastic  anode  foil  beyond  the  inner  cathode  radius  with  a  thin  metal  foil  (on  Gam¬ 
ble  11,  first  tantalum  then  aluminum  was  used11).  This  concept  has  been  tested  at  NRL  with  very 
favorable  results12  which  will  be  discussed  later.  The  numerical  simulations  reported  herein  modeled 
the  steady-state  operation  of  those  experimental  diodes.  The  empirical  data  showed  that  suppression  of 
the  peripheral  ion  emission  both  decoupled  the  diode  impedance  from  R/D  (thus  allowing  a  smaller  D 
for  a  given  impedance)  and  also  significantly  enhanced  the  ion  efficiency  in  comparison  to  a  diode  with 
the  same  impedance  but  with  larger  D  and  full  anode  ion  emission.  The  computer  runs  confirmed  the 
impedance  decoupling  from  R/D  but  failed  to  detect  any* Significant  ion  efficiency  enhancement.  How¬ 
ever,  the  numerically  observed  ion  efficiency  of  the  diode  with  ion  suppression  and  small  D  does  agree 
with  the  experimental  results.  Furthermore,  the  simulation  suggests  a  possible  physical  mechanism 
whereby  the  true  magnitude  of  the  net  ion  current  could  have  been  partially  masked  in  the  "bench¬ 
mark”  (i.e.,  full  anode  ion  emission)  experiments  by  a  comoving  electron  current. 

II.  THE  DIODE  AND  ITS  NUMERICAL  MODEL 

The  actual  Gamble  II  experimental  diode  configuration  is  drawn  in  Figure  2.  For  the  discussions 
presented  later  in  this  report,  it  is  particularly  important  to  note  the  structure  of  the  inner  cathode 
assembly.  The  "target"  is  a  sheet  of  KIMFOL  stretched  over  an  inner  cylindrical  cathode  shank  of 
radius  2.25  cm.  Only  the  ions  flowing  within  that  radius  can  be  extracted  from  the  diode  and  counted 
in  the  experimental  (/net^ok  measurements.  For  ease  of  numerical  modeling,  the  actual  diode  of  Fig¬ 
ure  2  was  simplified  slightly  to  the  form  depicted  in  Figure  3.  The  anode  foil  radius  was  shortened 
from  about  6.0  cm  to  3.7  cm.  Since  that  is  still  almost  a  centimeter  larger  than  the  cathode  outer  radius 
the  electric  field  structure  along  the  electron-emitting  cathode  shank  and  along  the  anode's  ion-emitting 
regions  opposite  the  cathode  shank  should  not  be  significantly  modified.  In  addition,  the  axially  bound¬ 
ing,  solid  cathode  and  anode  surfaces  are  taken  as  completely  planar  and  parallel.  The  front  surface  of 
the  cathode  "door"  (see  Figure  2)  coplanar  with  the  "target"  face.  Finally,  the  thickness  of  the  cathode 
shank  is  taken  as  0.35  cm  while  in  the  experiments  it  was  usually  0.25  cm.  The  extra  millimeter  allows 
for  somewhat  better  numerical  resolution  over  the  all-important  tip  of  the  shank  where  most  of  the  net 
electron  current  is  emitted.  Having  fixed  the  essential  geometry,  it  is  a  straightforward  process  to 
implement  the  numerical  model. 

NRL's  D10DE2D  computer  code  was  employed  to  numerically  simulate  the  steady-state  operating 
conditions  for  this  diode  for  various  sets  of  parameters.  The  details  of  the  code  may  be  found  else¬ 
where.13  It  is  sufficient  here  to  point  out  that  D10DE2D  calculates  equilibrium  electric  and  magnetic 
field  strengths  over  an  NZ  x  NR  mesh  of  discrete  data  points  on  a  predetermined  computational  region 
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Fig.  1  —  The  pinch-reflex  diode 


corresponding  to  an  arbitrary  R  —  Z  planar  cross-section  passing  through  the  diode's  centerline.  Com¬ 
plete  azimuthal  symmetry  is  assumed.  A  finite  number  of  macroelectrons  and  macroprotons  having 
correct,  physical  charge-to-mass  ratios  are  advanced  timestep-by-timestep  across  the  mesh  in  a  rela- 
tivistically  covariant  manner.  A  self-consistent,  steady-state  solution  is  sought  both  for  field  structures 
as  well  as  for  particle  flows.  No  time-dependent  phenomena  are  actually  treated. 

The  physical  volume  simulated  by  the  D10DE2D  code  extends  radially  outward  from  the  central 
axis  to  a  distance  of  6. IS  cm  in  the  vacuum  gap  and  axially  forward  from  the  plane  corresponding  to  the 
solid  cathode  face  out  to  the  opposing  plane  surface  of  the  solid  anode.  Since  it  presumes  azimuthal 
symmetry,  the  computer  code  only  deals  with  a  single  R  —  Z  planar  cross-section  extending  out  from 
the  centerline.  This  computational  region  is  presented  in  Figure  4.  The  grid  points  are  shown  as  dots 
and  correspond  to  the  centers  of  their  respective  rectangular  data  cells.  Counting  the  monolayer  of 
guard  cells  which  completely  surround  the  numerical  region  in  which  "particles"  are  "allowed,"  a  total  of 
(NZ  +  2)  x  (NR  +  2)  =  66  x  125  =  8,250  data  cells  are  used.  The  bottom  boundary  corresponds  to 
the  central  axis  ( R  =  0)  of  the  diode.  The  right  boundary  corresponds  to  the  solid  anode  plane  and  is 
maintained  at  the  full  anode  potential,  V.  The  left  boundary  is  kept  at  zero  voltage  since  it  represents 
the  solid  cathode  plane.  Finally,  the  upper  boundary  is  simply  a  source-free,  vacuum  interface  between 
what  may  be  loosely  described  as  the  "diode  region"  and  the  pulsed  power  transmission  line.  The  elec¬ 
trostatic  potential  is  therefore  graded  linearly  from  zero  to  V  along  that  upper  boundary.  Perfectly  con¬ 
ducting  cathode  and  anode  surfaces  are  included  in  the  computational  region  as  shown  in  Figure  4. 
Those  surfaces  are  treated  numerically  via  a  capacitance  matrix  technique  described  elsewhere.14  On 
the  order  of  I04  macroelectrons  and  104  macroprotons  participate  in  the  simulation  at  steady  state. 
These  macroparticles  are  emitted  at  their  respective  electrodes  along  the  heavy-lined  surfaces.  Axial 
currents  in  the  cathode  shank  and  in  the  anode  stalk  are  treated  rigorously  as  a  function  of  z  in  order  to 
ensure  an  accurate  distribution  of  Be  throughout  the  diode.  The  results  of  these  numerical  simulations 
as  well  as  a  summary  of  the  conclusions  which  can  be  drawn  from  them  are  presented  in  the  following 
two  sections. 

III.  SIMULATION  RESULTS 

In  order  to  test  the  theorized  effects  of  selective  ion  emission  suppression  on  the  performance 
characteristics  of  the  diode  of  Figure  2,  three  specific  test  cases  were  simulated.  In  all  three  cases,  the 
diode  voltage  was  fixed  at  1.8  MV  in  accordance  with  the  capabilities  of  the  Gamble  II  machine.  Cases 
I  and  II  were  both  run  with  an  A  -K  (anode-cathode)  gap  of  5  mm  just  as  shown  in  Figure  3.  Case  I 
corresponded  to  a  "benchmark"  run  in  which  no  suppression  of  ion  emission  on  the  anode  foil  was 
attempted.  Protons  were  emitted  along  the  anode  mesh  from  cells  UZ./R)  =  (48,  2)  through  (48,  65) 
just  as  indicated  by  the  heavy  line  in  Figure  4.  In  Case  II,  the  imposition  of  a  metal  foil  annulus  over 
the  outer  anode  face  was  simulated  by  shutting  off  ion  emission  beyond  cell  (48,  49).  Finally,  for  Case 
III,  the  restricted  ion  emission  was  retained  while  all  axial  dimensions  in  the  simulation  were  reduced 
by  a  factor  of  0.70.  This  changed  the  A  -K  gap  to  3.S  mm  and  reduced  the  diode  impedance  back  to  its 
value  in  Case  I.  The  detailed  results  of  these  simulations  are  as  follows. 

Case  I  was  run  4400  timesteps,  At,  of  1.5  x  10-12  seconds  each  until  an  equilibrium  state  was 
reached.  In  order  to  determine  the  true  steady-state  values  of  the  quantities  of  interest  a  new  statistical 
technique  was  applied  to  the  numerical  data.  Specie  current  values  were  tabulated  at  intervals  of  25A t 
over  the  last  200  steps.  These  values  were  then  subjected  to  the  same  type  of  least  squares  error 
analysis  and  data  reduction  methods  as  is  normally  reserved  for  the  treatment  of  experimental  observa¬ 
tions.  The  details  of  this  new  computational  approach  are  presented  in  the  Appendix.  The  equilibrium 
ion  and  electron  currents  for  this  "benchmark"  case  were  230  ±  6  kA  and  147  ±  6  kA  respectively. 
Given  the  diode  voltage  of  1.8  MV,  this  yields  an  impedance  of  4.77  ±0.11  ohms  for  the  device,  well 
within  the  desired  operating  regime  of  4-5  O.  The  net  ion  efficiency,  tj,  =  /«>„// diode  — 
W</io„  +  ^electron)'  is  then  0.390  ±  0.016.  This  is  in  excellent  agreement  with  the  40%  ion  efficiency 
recently  reported15  for  a  numerical  simulation  of  a  2.0  MV,  4  ohm  PRD  (pinch-reflex  diode)  having  an 
R/D  of  7.58  compared  with  an  R/D  of  6.0  for  this  run.  It  differs  significantly,  however  from  the  0.24 
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Fig.  4  _  Set-up  of  the  numerical  simulation  region 
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average  ion  efficiency  reported  from  a  series  of  Gamble  II  experimental  runs16  (see  Figure  5).  At  least 
part  of  this  discrepancy  between  the  numerical  and  experimental  results  can  be  explained  by  examining 
Figure  2.  The  location  of  the  /NET  Rogowski  coil  loop  in  the  experimental  set-up  shows  that  it  is 
measuring  only  that  portion  of  the  ion  current  which  hits  the  "target"  inside  the  cathode  shank.  Any 
ions  hitting  the  cathode  beyond  a  radius  of  2.25  cm  are  not  counted.  To  properly  compare  the  numeri¬ 
cal  results  to  the  experiment,  this  "effective”  ion  current  must  be  calculated.  Figures  6  and  7  present 
radial  profiles  of  the  emitted  and  collected  ion  current  densities  flowing  in  the  steady  state  diode. 
Integrated  from  these  profiles  are  not  only  the  total  ion  currents  but  also  the  effective  currents  falling 
within  R  =  2.25  cm  (indicated  by  the  dotted  lines).  This  shows  that  the  empirically  measured  number 
should  be  119  kA,  not  147  kA.  This  still  implies  an  ion  efficiency  of  0.316.  Although  much  closer  to 
the  experimental  data,  it  is  still  far  from  "good"  agreement.  A  more  complete  explanation  for  the  poor 
correlation  can  be  found  by  examining  the  nature  of  the  electron  currents  in  the  diode.  That  new 
phenomenon  will  be  discussed  after  the  numerical  results  for  Cases  II  and  III  are  presented. 

Having  completed  the  benchmark  run,  the  next  step  was  to  test  the  effect  of  restricted  ion  emis¬ 
sion  while  holding  all  other  diode  parameters  constant.  The  maximum  radius  for  ion  emission  along 
the  anode  foil  was  reduced  from  R  =  3.2  cm  to  R  =  2.4  cm.  Just  as  for  Case  I,  Case  II  was  run  to 
equilibrium  at  t  =  4400  timesteps.  At  that  point,  using  the  same  statistical  data  reduction  technique, 
the  currents  were  found  to  be  /,  =  186  ±  3  kA  and  /,  =  83  ±  2  kA.  These  figures  imply  a  diode 
impedance  of  6.69  ±  0.08  ohms  and  a  net  ion  efficiency  of  0.309  ±  0.007.  As  expected,  the 
impedance  has  been  very  significantly  raised  in  spite  of  the  constant  A  -K  gap  size,  D  =  5  mm.  Once 
again,  to  correlate  the  ion  efficiency  to  the  experimental  results,  the  radial  profiles  of  the  emitted  and 
collected  ion  current  densities  are  presented  in  Figures  8  and  9.  An  effective  /,  of  79  kA  is  then 
observed  to  strike  within  the  "target"  radius  of  2.25  cm.  Interestingly  enough,  the  effective  ion 
efficiency  has  dropped  a  mere  two  percentage  points  to  (r),)elr  =  0.294  in  spite  of  the  2  ohm  impedance 
rise. 


Finally,  for  Case  III  the  restricted  ion  emission  radius  of  2.4  cm  was  retained  while  the  gap  was 
reduced  from  5.0  to  3.5  mm.  The  observed  steady-state  diode  currents  were  /e  =  247  ±  3  kA  and 
/,  =  130  ±  4  kA  yielding  an  impedance  of  4.77  ±  0.07  ohms  and  a  total  ion  efficiency  of 
0.345  ±  0.012.  Thus,  as  expected,  the  impedance  of  Case  I  has  been  restored  in  spite  of  the  decreased 
gap  size,  D.  Although  it  had  been  hoped  that  the  ion  efficiency  would  increase,  a  drop  of  almost  0.05 
was  instead  recorded.  Hopes  that  a  comparison  of  the  effective  ion  currents  would  be  more  optimistic 
were  likewise  disappointed.  Figures  10  and  11  record  the  corresponding  Case  III  ion  emission  and  col¬ 
lection  profiles.  They  show  that  the  effective  ion  current  is  119  kA  —  exactly  the  same  as  for  Case  I. 
Thus  the  simulations  of  Case  I  and  Case  III  share  identical  impedances  and  identical  efficiencies.  This 
is  in  sharp  contrast  with  the  experimental  data.  The  same  diode  as  modeled  in  Case  III  fired  28 
separate  times  on  Gamble  II  yielded  the  modified  efficiency  data  shown  in  Figure  12  (overlapped  with 
the  Case  I  data).  With  an  average  effective  ion  efficiency  of  0.34,  it  represents  a  clear  departure  from 
the  experimental  data  garnered  from  toe  diode  without  emission  restriction  and  with  D  =  5.0  mm 
although  the  average  impedances  was  nearly  identical  for  the  two  cases.  On  the  other  hand,  the  Case 
III  efficiency  figure  of  0.34  is  in  reasonable  agreement  with  the  0.316  measured  numerically  for  both 
cases. 


As  previously  mentioned,  at  least  a  partial  explanation  for  the  disagreement  between  simulation 
and  experiment  may  rest  in  the  diode’s  electron  flow  patterns.  The  DIODE2D  computer  code  is  not 
electromagnetic;  rather,  it  is  electrostatic-magnetostatic.  Inductive  electric  fields  are  not  generated  by 
the  particle  currents  as  they  should  be  due  to  temporal  variations  in  B„.  This  limitation  allowed  sus¬ 
tained  electron  currents  of  huge  magnitude  to  arise  between  the  inside  cylindrical  surface  of  the  cathode 
shank  and  the  "target"  plane  inside  the  shank.  The  ability  of  the  ion  flow  in  a  diode  to  draw  in  elec¬ 
trons  to  neutralize  its  space  charge  is  well  established.'7  How  this  neutralizing  fill  of  electron  charge 
maintains  itself  in  the  A  -K  gap  is  not  well  understood.  At  least  in  this  simulation,  part  of  that  charge  is 
constantly  replenished  by  the  cathode-to-cathode  electron  currents  herein  observed.  For  Case  1,  this 
electron  current  comoving  with  the  ions  hitting  the  "target"  region  has  a  magnitude  of  almost  200  kA. 
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For  Case  (I  it  is  only  about  120  kA.  While  for  Case  111  it  increases  slightly  again  to  around  140  kA.  It 
is  probable  that  the  inclusion  of  a  physical,  inductive  f-fteld  would  greatly  reduce  the  size  of  those 
reverse  currents.  It  is  questionable  if  they  would  be  eliminated  entirely.  If  only  20  kA  of  electron 
current  were  left  comoving  whh  the  ion  beam  in  the  detector  region  for  the  Case  I  numerical  simula¬ 
tion,  the  experimentally  recorded  efficiency  of  0.24  would  be  achieved.  It  thus  seems  quite  possible 
that  such  comoving  electrons  may  be  masking  some  of  the  true  ion  current  being  extracted  from  the 
diode  in  the  Case  I  Gamble  II  experiment.  This  theory  remains  to  be  tested.  Figures  13,  14,  and  15 
provide  plots  of  sample  electron  positions  (each  particle  carries  the  same  predetermined  amount  of 
charge)  for  the  equilibria  of  Cases  I,  II,  and  III  respectively.  They  show  how  electron  space  charge  fills 
the  entire  interior  of  the  diode.  Unfortunately,  they  do  not  show  explicit  flow  patterns  inside  the 
cathode  radius.  One  thing  that  does  seem  clear  is  the  difference  in  cathode-tip-to-anode-foil  electron 
flow.  For  Case  I,  without  ion  emission  suppression,  the  ion  space  charge  in  the  gap  between  the  tip  and 
the  foil  allows  electrons  to  stream  almost  directly  across,  thus  illuminating  the  foil  face  out  to  about  3 
cm  with  significant  negative  space  charge.  Electrons  hugging  that  surface  out  to  large  radii  encourage 
enhanced  ion  emission.  On  the  other  hand,  Figures  14  and  15  show  that  elimination  of  the  ion  space 
charge  fill  near  the  cathode  tip  causes  much  of  the  electron  stream  there  to  pinch  radially  inward  early 
in  its  transit  of  the  gap.  Thus  the  mainstream  of  the  electron  flow  on  the  average  for  Cases  II  and  III 
stays  further  away  from  the  foil  face  than  in  Case  I.  This  reduces  expected  ion  emission  and  can 
apparently  reduce  the  ion  efficiency  benefits  of  decreasing  the  gap  size,  D. 

IV.  CONCLUSIONS 

From  the  results  summarized  in  Table  1,  one  may  conclude  that  selective  ion  emission  suppres¬ 
sion  definitely  allows  the  use  of  smaller  A  -K  gaps  (therefore,  higher  R/D  ratios)  for  a  fixed  diode 
impedance.  This  can  be  beneficial  for  all  light  ion  beam  work  in  which  efficient  generator-to-diode 
power  transmission  dictates  a  high  diode  impedance  since  a  smaller  A-K  gap  means  that  the  ions  must 
spend  less  time  subject  to  the  trajectory-bending  effects  of  the  strong  electric  and  magnetic  fields  found 
in  the  diode.  The  resultant  impact  on  ion  beam  quality  is  illustrated  in  Figures  16,  17,  and  18  which 
plot  sample  proton  equilibrium  positions  for  Cases  I,  II,  and  III  respectively.  The  ions  are  emitted  only 
at  data  cell  centers  along  the  anode  face.  Thus  they  are  initially  spaced  out  evenly  in  coaxial,  laminar 
streams.  The  relative  smoothness  of  each  stream  may  be  followed  across  the  gap  by  tracing  from  point 
to  point.  Significantly  less  erratic  stream  behavior  is  observed  in  Figure  18  (Case  III)  with  its  3.5  mm 
gap  than  for  either  of  Cases  I  and  II  in  Figures  16  and  17.  Although  this  confirms  at  least  one  positive 
effect  of  partial  ion  emission  suppression,  the  simulations  did  not  show  any  decoupling  of  the  effective 
ion  production  efficiency  from  the  diode  impedance  as  was  found  in  the  corresponding  Gamble  II 
experiments.  Although  the  simulation  results  are  clouded  by  unrealistically  large  electron  currents 
comoving  with  the  ion  beam,  the  physical  existence  of  electron  currents  even  an  order-of-magnitude 
smaller  than  those  numerically  observed  would  be  sufficient  to  explain  the  discrepancy  between  experi¬ 
ment  and  simulation. 

Two  distinct  elements  of  follow-on  research  should  be  undertaken.  First  of  all,  the  simulations 
described  herein  should  be  repeated  using  a  fully  electromagnetic  computer  code.  The  PREMAS  com¬ 
puter  code  (the  successor  of  DIODE2D  at  NRL)  is  currently  being  upgraded  to  include  time  dependent 
electric  and  magnetic  field  solving  algorithms.  It  can  be  applied  to  this  problem  as  soon  as  that  coding 
and  debugging  is  complete.  Secondly,  an  experimental  diagnostic  technique  to  separate  out  the  pure 
ion  beam  current  from  any  comoving  electron  component  should  be  designed.  The  Gamble  II  runs 
could  then  be  repeated  with  such  a  new  diagnostic  in  place. 
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radial  profile  of  collected  ion  current  density  Tor  Case  I 


Fig.  8  —  Radial  profile  of  emitted  ion  current  density  for  Case  II 
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Fig.  10  —  Radial  profile  of  emitted  ion  current  density  for  Case  III 
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Fig.  1 1  —  Radial  profile  of  collected  ton  current  density  for  Case  III 
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Fig.  15  —  Sample  steady-state  electron  plot  for  Case  III 
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To  illustrate  the  data  reduction  procedure,  the  ion  current  of  Case  II  in  this  report  will  be 
analyzed.  The  last  eight  collected  ion  current  values  (25  At  apart)  were  found  to  be  83,  81,  87,  84,  82, 
78,  93,  and  80  kA  respectively.  This  yields  a  mean  value,  <  /,COii>  ,  of  83.5  kA.  To  determine  the  ran¬ 
dom  sampling  error  associated  with  that  value,  the  mean  value  of  the  square,  <  /,2coii  ^  >  ,s  found  to  be 
6991.5  kA2  from  which  the  variance, 

O’icoll  “  >/  <  (/coll  >  —  <  A  coll  >  2  (A.l) 

is  computed  to  be  4.387  kA.  The  standard  deviation  of  the  mean  of  the  given  eight  values  is  then  sim¬ 
ply 

(®’/coii)<n  ”  o' (coil  “  1.55  kA.  (A. 2) 

Thus,  on  the  basis  of  the  computer-supplied  data,  the  collected  ion  current  for  Case  II  can  be  written  as 
( 7() coii  ”  83.5  ±  1.6  kA.  In  exactly  the  same  manner,  the  last  eight  emitted  ion  current  values  of  97, 
99,  112,  73,  48,  85,  77,  and  76  kA  reduce  to  a  steady-state  prediction  of  “  83.4  ±  6.5  kA. 

The  much  larger  standard  deviation  stems  directly  from  the  higher  residuals  of  each  data  point  with 
respect  to  their  mean  value.  To  arrive  at  a  specific  prediction  for  the  net  equilibrium  ion  current 
flowing  through  the  diode,  the  weighted  average  is  taken  of  the  emitted  and  collected  ion  currents  using 
the  method  of  least  squares  with 

.  ^  Amit^/A^emil)*?  ^  Aaall^*/ (<rco||)m  , 

>«  +‘l/(o-Co,|U 


*1  (<remil)2  T  ' 

Applying  these  formulae  lead  to  the  result  quoted  in  the  text  of  this  report,  /jon  -  83  ±  2  kA. 

An  identical  analysis  is  then  applied  to  find  the  electron  current  and  its  error  brackets.  The  net 
diode  current  is  then  simply 

Aliode  “  Aon  4"  ^electron  (A. 5) 


(tf-diode)*  ”  (O'iofi)m  +  (^electron) l-  (A.6) 

From  this  quantity  flows  the  expressions  for  the  diode  impedance,  z,  and  ion  production  efficiency 
For  the  impedance,  one  may  write  z  —  K//djode  where  V  is  the  constant  diode  voltage.  The  standard 
deviation  from  the  mean  of  this  quantity  is  then  given  by 


<*.>*“  7T 


(''’diode)  m- 


Similarly,  the  ion  production  efficiency,  t\,  -  Aon/ Anode*  has  a  standard  deviation  found  from  the  expres¬ 
sion, 

+  i^r.  <a.») 

I  'diode  I  'diode 


(O'diode)m  + 


The  above  equations  summarize  the  error  analysis  which  flows  naturally  from  the  reduction  of  raw, 
nutnerical  data.  It  provides  a  simple  framework  for  the  statistically  valid  presentation  of  steady-state 
computer  simulation  code  results  for  diode  physics  research. 


STATISTICAL  TREATMENT  OF  COMPUTATIONAL  DATA 


Whenever  numerical  predictions  extracted  from  computer  simulations  are  cited  in  the  literature 
they  are  almost  always  presented  without  any  indication  of  the  error  brackets  associated  with  those 
predictions.  Conventional  error  analysis  such  as  that  performed  routinely  for  experimental  data  is  tradi¬ 
tionally  ignored  by  many  computational  physicists.  When  one  analyzes  a  given  physical  system  using  a 
computer  simulation  code,  the  numerical  data  distilled  from  the  code  is  just  as  prone  to  systematic  and 
random  sampling  errors  as  that  gleaned  from  a  physical  experiment.  Programming  flaws  in  the  code  as 
well  as  physical  flaws  and  resolution  limits  in  the  simulation  model  will  introduce  systematic  errors  into 
the  resulting  data.  Similarly,  when  dealing  with  a  computer  code  which  attempts  to  find  the  steady-state 
operating  conditions  of  a  given  system,  periodic  sampling  of  the  parameters  of  interest  as  they  relax 
into  their  equilibrium  values  will  give  rise  to  random  errors  around  the  "true*  values.  Both  types  of 
errors  can  be  quantified  in  order  to  allow  an  outside  observer  to  know  how  much  weight  he  should  give 
to  the  computational  results. 

The  identification  and  minimization  of  gross  systematic  errors  contained  in  the  algorithms  of  a 
given  computer  code  is  normally  accomplished  once  and  for  all  before  a  given  code  is  accepted  for  use 
in  production  simulation  runs.  Careful  debugging  using  a  variety  of  standard  test  cases  can  eliminate 
numerical  deviations  from  true  physics  that  go  beyond  established  assumptions  for  the  systems  to  be 
modeled.  In  addition  to  computer  programming  bugs,  there  are  also  the  finite  resolution  limitations 
inherent  in  any  digitized  model.  Some  of  these,  such  as  standard  numerical  "round-off"  errors  can  be 
made  insignificant  through  clever  coding  of  the  numerical  algorithms.  Other  sources  of  error  such  as 
the  finite  number  of  spatial  data  points  and  the  finite  number  of  particles  used,  give  rise  to  actual 
"minimum  errors"  which  can  have  significant  sizes  compared  with  the  values  of  interest  and  which  can 
only  be  lowered  by  changing  the  numerical  parameters  or  by  rewriting  sections  of  the  computer  code 
itself.  Any  potential  user  of  the  simulation  results  must  rely  upon  the  professional  abilities  of  the  indi¬ 
vidual  computational  physicist  regarding  the  removal  of  all  "bugs"  from  his  code.  On  the  other  hand, 
the  user  must  be  aware  that  there  will  definitely  be  some  systematic  errors  due  to  numerical  discrete¬ 
ness. 


The  treatment  of  random  sampling  errors  is  less  a  question  of  computational  talents  then  of 
scientific  rigor.  Observationalty  obtained  data  is  inherently  bracketed  by  error  bars.  Only  a  finite 
number  of  samples  of  data  can  be  obtained  from  a  given  system.  Like  its  physical  counterpart,  a 
numerical  system  cannot  be  truly  static.  It  will  oscillate  around  an  equilibrium  state  and  it  is  at  various 
points  on  these  "ripples"  that  the  diagnostic  data  is  output  from  the  program.  Therefore,  random  sam¬ 
pling  errors  must  occur  and  they  can  be  quantified  using  conventional,  statistical  treatments.  The  data 
reduced  for  presentation  in  this  report  may  be  used  to  illustrate  the  procedure.  The  major  macroscopic 
quantity  monitored  in  these  simulations  is  the  net  current  flowing  through  the  diode.  This  current  is 
made  up  of  an  electron  and  an  ion  component.  For  each  specie  component,  the  emitted  current  and 
the  collected  current  are  measured  separately.  Therefore,  there  are  four  "subcurrents"  that  are  directly 
monitored  and  must  be  analyzed  separately.  In  a  given  DIODE2D  simulation  run,  the  subcurrents  and 
net  charges  in  the  system  are  tallied  every  25  timesteps.  After  some  specific  "risetime”  is  passed,  these 
quantities  begin  to  oscillate  around  certain  plateau  values.  ■  Such  behavior  signals  the  system's  arrival  at 
an  equilibrium  state.  AJ  that  point,  the  last  eight  subcurrent  data  values  spanning  the  last  200 
timesteps  are  reduced  to  arrive  at  approximate  steady  state  predictions  for  the  diode  being  modeled. 


For  this  specific  problem  of  light  ion  diode  error  analysis,  another  element  of  statistical  analysis, 
namely  —  curve  fitting,  was  employed.  Its  purpose  was  to  obtain  the  best  possible  analytic  approxima¬ 
tion  for  the  radial  profile  of  the  ion  current  density  exiting  the  diode.  Such  an  explicit  representation 
for  7(r )  is  useful  particularly  for  theoreticians  attempting  to  model  the  transport  of  the  emerging  ion 
beam  away  from  the  diode  through  a  predetermined  transport  channel.1819  Previous  simulations  have 
observed  an  approximate  inverse  proportionality  of  the  current  density  to  the  radius.20  A  functional 
form  of 

7(r)  =  J0  +  air  -  r0)_1  (A.9) 

is  therefore  chosen  as  the  starting  point.  The  raw  data  of  averaged  ion  emission  and  collection  at  each 
mesh  point  in  the  r-dimension  along  the  anode  and  cathode  respectively  are  first  plotted  as  single  points 
(see  Figures  6-11).  Then  a  curve  of  form  (A.9)  is  fitted  to  it.  In  order  to  obtain  that  curve,  the 
method  of  least  squares21  is  employed.  The  three  unknowns  in  the  equation  are  then  given  by 


£  J,r2  T  r,  £  J,r, 

i  j,n  n  i  y, 

I  J}rt  1  J-  I  7, 2 


(A.  10) 


where 


Lr(2  I  r,  I  Jj2 
r0  -  D1  Lr,  NT  J,r, 

T  in  s  1  t  J}n 

t  n2  £  Jj2  t  J,n 

■  J0r0  +  d  '  i  n  £  in  i  J, 

T  Jjfj  £  Jjf,  T  7 ? 

I  n1  i  n  i  j,r, 

D  -  T  n  N  TJ, 

T  J,n  1  J<  tj? 


and  where  the  summations  are  from  /  »  1  to  i  —  N  (N  being  the  number  of  radial  data  points),  r,  is  the 
radius  of  the  "i"th  point,  and  Jt  is  the  observed  7(/j).  The  "goodness"  of  the  fit  of  the  data  to  the 
curve  thus  generated  can  be  gauged  by  comparing  the  root  mean  square  deviation  of  the  data  points 
from  the  curve  to  the  root  mean  square  deviation  of  those  same  points  from  their  average  value 
<7>  -  N~l  T  Jt.  The  resultant  "correlation  coefficient’  may  then  be  expressed  as22 


T  U  -  70-  a  (n~  r0)-') 
T  (Jt  -  <J» 2 


(All) 


It  should  be  noted  that  in  practice  some  care  must  be  exercised  in  the  application  of  (A.  11)  to  the 
problem  in  question.  Near  r  —  0,  the  values  of  7  become  quite  large  so  that  deviations  from  the 
smoothed  fit  tend  to  be  exaggerated.  To  avoid  that  distortion,  the  data  at  r  -  rx  is  ignored  in  the  calcu¬ 
lation  of  (A.l  1)  as  presented  in  Figures  6  through  11. 
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